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TAUOLA for simulation of tau decay and production: perspectives for 
precision low energy and LHC applications 

Z. Was a * 

a Institutc of Nuclear Physics, Polish Academy of Sciences, 
ul. Radzikowskiego 152, 31-342 Cracow, Poland 

The status of Monte Carlo system for the simulation of r-lepton production and decay in high-energy accelerator 
experiments is reviewed. Since previous r-lepton conference in 2008 some practical modifications have been 
introduced: (i) For the TAUOLA Monte Carlo generator of r-lepton decays, automated and simultaneous use of 
many versions of form-factors for the calculation of optional weights for fits was developped and checked to work 
in Belle and BaBar software environment. Work on alternative paramterizations of hadronic decays is advanced, 
(ii) the TAUOLA universal interface based on HepMC (the C++ event record) is now public. A similar interface 
for PHOTOS is now also public, (iii) Extension of PHOTOS Monte Carlo for QED bremsstrahlung in decays featuring 
kernels based on complete first order matrix element are gradually becoming widely available thanks to properites 
of the new, HepMC based interface, (iv) Tests of the programs systematized with the help of MC-TESTER are now 
available for FORTRAN and C++ users. 

Presented here results illustrate the status of the projects performed in collaboration with Nadia Davidson, 
Piotr Golonka, Gizo Nanava, Tomasz Przedzihski, Olga Shekhovtsova, Elzbieta Richter-Was, Pablo Roig, Qingjun 
Xu and others. 
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1. Introduction 



The TAUOLA package [1121314] for the simula- 
tion of r-lepton decays and PHOTOS [51617] for the 
simulation of QED radiative corrections in de- 
cays, are computing projects with a rather long 
history. Written and maintained by well-defined 
(main) authors, they nonetheless migrated into 
a wide range of applications where they became 
ingredients of complicated simulation chains. As 
a consequence, a large number of different ver- 
sions are presently in use. Those modifications, 
especially in case of TAUOLA, are valuable from 
the physics point of view, often did not found the 
place in the distributed versions of the program. 
Even if from the algorithmic point of view, ver- 
sions differ only in details, they incorporate many 
specific results from distinct r-lepton measure- 
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ments or phenomenological projects. Such ver- 
sions were mainly maintained (and will remain 
so) by the experiments taking precision data on 
r leptons. Interesting from the physics point 
of view changes are expected to be developped 
in FORTRAN. That is why, significant part of the 
TAUOLA should remain in FORTRAN for the forth- 
coming several years. 

However many new applications were devel- 
oped recently, often requiring a program inter- 
face to other packages (e.g. generating events for 
LHC, LC, Belle or BaBar physics processes). For- 
tunately, co-existence with CH — h is not a problem, 
at least not from the software point of view. 

The programs structure, prepared for the con- 
venience of FORTRAN users, was presented during 
previous r conferences, and we will not repeat it 
here. This time, let us concentrate on new in- 
terfaces for applications based on HepMC [8\ event 
record. We will also report on progress in im- 
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plementation (re-implementation) of techniques 
useful for fits. Analyses of high precision, high 
statistic data from Belle and BaBar are expected 
to progress from these solutions. 

Our presentation is organized as follows: Sec- 
tion 2 is devoted to new interfaces of TAUOLA and 
PHOTOS to HepMC applications of C++, where gen- 
uine weak and transverse spin effects can be taken 
into account as well. Section 3 is devoted to the 
discussion of optional weights in TAUOLA and their 
use for fits to data at the level of comparison with 
raw data. Progress on the work on new currents 
for hardonic decays which can be confronted with 
(tuned to) data using such optional weights en- 
abling simultaneous control of all experimental 
effects is mentioned too. In section 4 we present 
some new results for PHOTOS Monte Carlo for ra- 
diative corrections in decays. Section 5 is devoted 
to MC-TESTER, the program which can be used for 
semi-automatic comparisons of simulation sam- 
ples originating from different programs. 

Because of the limited space of the contribu- 
tion, and sizable amount of other physically in- 
teresting results, some of them will be excluded 
from conference proceedings. They will find place 
in future works, possibly with collaborators men- 
tioned in the Abstract. For these works, the 
present paper may serve as an announcement. 

2. TAUOLA universal interface and PHOTOS 
interface in C-\ — \- 

In the development of packages such as TAUOLA 
or PHOTOS, questions of tests and appriopriate 
relations to users' applications are essential for 
their usefulness. In fact, user applications may 
be much larger in size and human effort, than the 
programs discussed here. Good example of such 
'user applications' are complete environments to 
simulate physics process and control detector re- 
sponse at the same time. Distributions of final 
state particles distributions are not always of di- 
rect interest. Often properties of intermediate 
states such as spin state of r-lepton, coupling 
constants or masses of intermedate heavy parti- 
cles are of prime interest. As a consequence it is 
useful that such intermediate state properties arc 
under direct control of the experimental user and 




Figure 1. Scheme of Monte Carlo simulation system 
with communication based on event record. 



can be manipulated to understand detector re- 
sponses. Our programs worked well with FORTRAN 
applications where HEPEVT event record was used. 
Now, for the C++ HepMC [5] case, interfaces had to 
be re-written, both for TAUOLA [TT] and PHOTOS 
[12] , The interfaces are gradually enriched and 
for example for the PHOTOS improved kernel, is 
already available for the general use, as explained 
in ref. [T3]. A complete (not longitudinal only) 
spin correlations are available for Z/"f* decay in 
TAUOLA interface now. They are based on elec- 
troweak corrections taken from refs. [14115] . See 
Fig. [T] for the scheme of programs communica- 
tions. The size of the electroweak effects on lon- 
gitudinal r polarization is shown in Fig. [5] Such 
modular organization opens ways for efficient al- 
gorithms to understand detector systematics, but 
at the same time responsability to control soft- 
ware precision must be shared by the user. For 
that purpose automatization of tests MC-TESTER 
was prepared [9]. New functionalities were intro- 
duced into the testing package [10 . In particular, 
it works now with the HepMC event record, the 
standard of C++ programs and the spectrum of 
available tests is enriched. For details see Section 
5. 
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Figure 2. Polarization of r leptons produced from 
up quarks (Fig. (a)) and down quarks (Fig. (b)) at 
cos 6 = —0.2 (qq — > t + t~). For the red line gen- 
uine weak corrections are taken into account, the black 
represent Standard Born (default in the interface). 
The blue line is Born according to alpha scheme of 
\14\15f . The small bump on the red line on Fig. (a) 
is due to the WW threshold. 



3. Optional weights for TAUOLA Monte 
Carlo 



On one hand high statistics experimental samples 
are collected, and the obtained precision is high, 
on the other hand, there is a significant cross- 
contamination between distinct r decay channels. 
Starting from a certain precision level all chan- 
nels need to be analyzed simultaneously. Change 
of parametrization for one channel contribution 
to the background may be important for the fit 
of another one. This situation leads to a com- 
plex configuation where a multitude of parame- 
ters (and models) needs to be simultaneously con- 
fronted with a multitude of observables. One has 
to keep in mind that the models used to obtain 
distributions in the fits may require refinements 
or even substantial rebuilds as a consequence of 
comparison with the data. The topic was covered 
in detail in the r Section of Ref. [16] earlier this 
year. 

From the statistical point of view it is best to 
resolve such system in one automated step us- 
ing for example a method such as |17ll8j . This 
can be of course very dangerous from the point 
of view of systematic error control. But we will 
not elaborate on this point any further. From 
the technical point it is necessary, to calculate for 
each generated event (for each present in it: decay 
of t + and t~ separately) alternative weights; the 
ratios of matrix element squared obtained with 
new currents provided by the user, and the one 
actually used in generation. Then the vector of 
weights can be obtained and used in fits. We 
have checked that such a solution not only can 
be easily installed into TAUOLA as a stand-alone 
generator but it can also be incorporated into the 
simulation frameworks of Belle and BaBar col- 
laborations rather easily. The scheme as given in 
Fig. [3] was shown to work in real conditions. 

At the present step of our work no new mod- 
els of hadronic currents are prepared. Work on 
currents based on refs. [19|20|21] is taking shape 
mainly thanks to efforts by Olga Shekhovtsova. 
A dedicated patch is being prepared, to be later 
installed into users software environments. 



Physics of t lcpton decays requires sophis- 
ticated strategies for the confrontation of phe- 
nomcnological models with experimental data. 
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Figure 3. Flow chart for fifo communication and 
reweigting TAUDLA matrix elements, verified to be 
compatible with Belle and BaBar software. 



4. PHOTOS Monte Carlo for bremsstrahlung 
and its systematic uncertainties 

Thanks to exponentation properties and fac- 
torization, the bulk of the final state QED 
bremsstrahlung can be described in a universal 
way. Howevers, the kinematical configurations 
caused by QED bremsstrahlung are affecting in 
an important way signal/background separation. 
It may affect selection criteria and background 
contaminations in quite complex and unexpected 
ways. In many applications, not only in r decays, 
such bremsstrahlung corrections are generated 
with the help of the PHOTOS Monte Carlo. That 
is why it is of importance to review the precision 
of this program as documented in Refs. |5|6|7j . 
For the C++ applications, the first version of the 
program is also available now, it is documented 
in Ref. Q2]. 

In CH — I- applications, the complete first order 
matrix element for the two body decays of the Z 
into pair of charged leptons [13] is now available. 
Kernels with complete matrix elements, for the 
decays of scalar B mesons into a pair of scalars 
[12], for W ± ->• ^7 and for 7* -> tt+tt" [2J 



are available for tests or specially oriented decay 
particles frames. It will be rather easy now to 
integrate NLO kernels into the main version of 
the program, because of better control of decay 
particle rest frame than in the FORTRAN interface. 

In all of these cases the universal kernel of 
PHOTOS is replaced with the one matching exact 
first order matrix element. In this way terms 
necessary for NLO/NLL precision level are im- 
plementecQ. A discussion relevant for control of 
program systematic uncetainty in t — > -kv decay 
can be found in Ref. |25j . 

The algorithm covers the full multiphoton 
phase-space and becomes exact in the soft limit. 
This is rather unusual for NLL compatible algo- 
rithms. One should not forget that PHOTOS gener- 
ates weight-one events too, and does not require 
any phase space ordering. There is a full phase 
space overlap between the one where hard ma- 
trix element is used and the one for iterative pho- 
ton emissions. All interference effects (between 
consecutive emissions and emissions from distinct 
charged lines) are implemented with the help of 
internal weights. 

The results of all tests of PHOTOS with a NLO 
kernel confirm sub-permille precision level. This 
is very encouraging, and points to the possible 
extension of the approach outside of QED (scalar 
QED). In particular to the domain of QCD or if 
phcnomenological lagrangians for interactions of 
photons need to be applied. For that work to 
be completed spin amplitudes need to be further 
studied. Let us point to Ref. [26] as an example. 

Kernels necessary for NLO are only available 
at present as options for tests only. They can be 
uploaded from the PHOTOS web page [27]. In the 
C++ applications they will be gradually installed 
into the main version of the program. The first 
step was the case of Z/j* decay and is already 
completed now. This opens the way to discuss 

1 Notc that here the LL (NLL) denote collinear logarithms 
(or in case of differential predictions terms integrating into 
such logarithms). The logarithms of soft singularities are 
taken into account to all orders. This is resulting from 
mechanisms of exlusive exponentiation |24| of QED. The 
algorithm used in PHOTOS Monte Carlo is compatible with 
exclusive exponentiation. Note that our LL/NLL precision 
level would even read as respectively NLL/NNNLL level 
in some naming conventions of QCD. 
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systematic errors in user denned acceptance re- 
gions. This will be an important supplement to 
tests-comparisons with other programs such as 
KKMC [2S], see Fig. gj 

5. MC-TESTER and user denned tests. 

Our work on MC-TESTER [10] reached maturity. 
As in the past, the program main purpose remain 
benchmarking the decay part of different Monte 
Carlo chains. Generated events have to be stored 
in event records: be it of FORTRAN or C++. De- 
fault distributions consist of all possible invariant 
masses which are automatically generated and 
stored for each found decay channel of the par- 
ticle under test. Then, at the analysis step, in- 
formation from a pair of such runs may be com- 
pared and represented in the form of tables and 
plots. At present, users macros can be easily in- 
stalled, in particular all demo distributions given 
in papers on C++ interfaces for TAUOLA [II] and 
PHOTOS [12] can be obtained in that way. Set-up's 
for benchmarking the interfaces, such as interface 
between r-lepton production and decay, including 
QED bremsstrahlung effects can be prepared in 
that way. 

The updated version of MC-TESTER was found 
useful for FORTRAN [7IT5] and for C++ fH] exam- 
ples where spurious information (on soft photons) 
was removed. 

Finally, let us mention that the program is 
available through the Grid Project LCG/Genser 
web page, see Ref. [29] for details. This is the case 
for TAUOLA C++ and for PHOTOS C++ in the near 
future as well. The FORTRAN predecessors have al- 
ready been available for some time. 

6. Summary and future possibilities 

The status of the computer programs for the 
decay of r leptons TAUOLA and associated projects 
TAUOLA universal interface and MC-TESTER 
was reviewed. The high-precision version of 
PHOTOS for radiative corrections in decays, was 
presented also. All these programs are available 
now for C++ applications thanks to the HepMC 
interfaces. 

New results for PHOTOS were mentioned. For 



the Z decay channel the complete next-to-leading 
collinear logarithms effects can now be simulated 
in C++ applications. However, in most cases 
these effects are not important, leaving the stan- 
dard version of PHOTOS sufficient. The exam- 
ple is important for the general case, it helps to 
better understand questions related to matching 
hard emission matrix elements with parton show- 
ers without the necessity to introduce any bound- 
aries within the phase space. Thanks to this work 
the path for fits to the data of electromagnetic 
form-factors is opened. The effects of these form- 
factors may be significantly larger and physically 
more interesting than complete next-to-leading 
order effects of QED or scalar QED. 

The presentation of the TAUOLA general- 
purpose interface in C++ was given. It is now 
more refined than the FORTRAN predecessor. Elec- 
troweak corrections can be used in calculation of 
complete spin correlations in Z/^* mediated pro- 
cesses. 

Distinct versions of the hadronic currents for 
TAUOLA library are available for the FORTRAN 
distribution. Work on new ones and on refine- 
ments of old strategies of fits using alternative 
weights and projection operators is progressing. 

The new version of MC-TESTER is stable now. 
It not only works with HepMC of C++ but enables 
user defined tests in experiments' software envi- 
ronments. 
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Figure 4. The spectrum of the 77 invariant mass 
in Z — > fi + fi^n^ decay. Events with two hard 
photons, both of energy above 1 GeV in the Z 
rest frame are taken and the invariant mass of the 
photon pair, normalized to Z mass is shown: for 
CEEX2 and CEEX1 (case a), and CEEX2 and 
PHOTOS (case b). The prediction from PHOTOS is 
clearly superior for applications aiming at simula- 
tions for Higgs boson backgrounds than CEEX1. 
In the case of solution based on YFS exponen- 
tiation \2ty . the second order matrix element 
must be taken into account. Fig. b was 06- 
tained with the help o/examples/testing/Zmumu 
PHOTOS demonstration example ( as documented 
in Ref. fWj). 



